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We investigate the Hagedorn transitions of string networks with Y-junctions as may occur, for 
example, with (p, q) cosmic superstrings. In a simplified model with three different types of string, 
the partition function reduces to three generalised coupled XY models. We calculate the phase 
diagram and show that, as the system is heated, the lightest strings first undergo the Hagedorn 
transition despite the junctions. There is then a second, higher, critical temperature above which 
infinite strings of all tensions, and junctions, exist. Conversely, on cooling to low temperatures, only 
the lightest strings remain, but they collapse into small loops. 



I. INTRODUCTION 

The statistical mechanics of string networks has been 
the object of numerous studies because of the importance 
of strings or string-like entities across all energy scales. 

In general, either because of the large number of con- 
figurational microstates or because of the large number 
of excited quantum states that such a network possesses, 
the networks undergo transitions in which, as tempera- 
tures rise, strings proliferate. In the language of config- 
urational states such a transition is termed a Feynman- 
Shockley transition, after Feynman's description of the A- 
transition of He in terms of vortex production pQ . From 
the viewpoint of counting excited states it is called a 
Hagedorn transition [2]. [Henceforth we follow the com- 
mon usage of Hagedorn transition to apply to both cases, 
which are similar in structure in many ways.] 

Specifically, in QCD, the sudden proliferation of colour 
flux tubes (the original dual hadronic strings) explains 
quark deconfincmcnt as temperature rises (see, for exam- 
ple, [31I3J|S]). I n cosmology at the GUT scale, where cos- 
mic strings arise in all reasonable supersymmetric mod- 
els incorporating electroweak unification [BJ, the statis- 
tical mechanics of cosmic string networks has been in- 
vestigated in order to understand their properties at for- 
mation and their late time scaling solutions, crucial for 
determining their cosmological consequences [7J [5] . For 
fundamental strings there has been substantial work on 
exploring the effects of such transitions on the extremely 
early universe [TOJ HU EE2J E3] . 

More recently, attention has turned again to funda- 
mental string networks, following new developments in 
superstring theory. Indeed, a network of cosmic super- 
strings is expected to form when a brane and anti-brane 
annihilate at the end of string-motivated brane inflation 
models. The network contains fundamental F-strings, 
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Dirichlct D-strings, and (p, g)-strings which are bound 
states of p F-strings and q D-strings [HI H21 [TBJ HZ], 
meeting at Y-junctions (or vertices). The presence of 
Y-junctions, as well as the spectrum of tensions of the 
strings, is a key characteristic of such networks and leads 
to more complicated dynamics. Much work has been 
done to determine how (j>, q)-like string networks evolve, 
both by analytic methods and numerical simulations, 
with particular regard to scaling solutions, their effect 
on the CMB as well as other observable consequences 
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Other than being stable against break-up, such strings 
differ from earlier superstrings in that, due to the warp- 
ing of space-time, their tensions are not of the Planck 
scale but many orders of magnitude smaller. As a result 
any Hagedorn transitions may even arise later than the 
reheating of the universe, and hence be of direct relevance 
for astrophysics. A necessary first step in seeing whether 
this is the case is to determine the phase diagram for the 
Hagedorn transitions of a network with more than one 
type of string, and this is the goal of the present paper. 

Our approach is to attempt to map the thermodynam- 
ics of string networks with junctions into the thermody- 
namics of a set of interacting dual fields, whereby the 
Hagedorn transitions of the strings become conventional 
transitions of the fields, a situation with which we are fa- 
miliar. One can imagine several ways to attempt this. We 
adopt the simplest, generalising the methods for describ- 
ing quark deconfincmcnt mentioned above (with its flux- 
tube Y-junctions) to something more like (p, g)-strings. 

Hence we investigate the equilibrium statistical me- 
chanics of cosmic superstring networks using methods 
motivated by [21 IH [5] . However, it is important to note 
that there is at least one major difference between cosmic 
superstrings and QCD fluxlines: with multiple tensions 
(from different string types), we expect cosmic super- 
string networks to show multiple Hagedorn transitions. 

In subsequent sections we derive and analyse the phase 
structure of a f/iree-string model with junctions. This 
is a reduced model of realistic cosmic superstrings, for 
which (p, q) G Z x Z form a doubly-infinite family. Since 
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FIG. 1: Different critical temperatures for our simplified 
model of cosmic superstrings (with tensions o\ < a 2 < 03) 
with Y-junctions. The lower Hagedorn temperature Ti is de- 
termined by <7i whereas the higher Hagedorn temperature T, 
is a determined by all the a a (a = 1,2,3). n v denotes the 
density of vertices (or Y-junctions) joining infinite strings at 
temperature T. 



string tension (or energy/unit length) increases with p, q, 
all but low values will be suppressed at high tempera- 
ture. We therefore adopt the simplest non-trivial scheme, 
taking the two lightest strings and their bound-state 
(and anti-strings), all which have different tensions <r a , 
a = (1,2,3). For example, depending on parameters, 
these could be the (1,0), (0,1) and (1,1) strings. We 
show that as the system is heated, the lightest tension 
strings first undergo the Hagedorn transition, despite the 
presence of Y-junctions. Conversely, at low tempera- 
tures, only the lightest strings remain, before they dis- 
appear into loops. Our results are summarized in figure 

m 

These conclusions may have important consequences 
for (p, q) string networks in that, if only the lightest 
strings remain after a non-adiabatic quench, no signif- 
icant role would be played by the junctions whose prop- 
erties have been studied so extensively. The dynamics 
would then be that of a single string type with no junc- 
tions (though there may be loops containing strings of 
different types; as explained below, our analysis is limited 
to infinite strings). This is not an idle proposition in that, 
although our analysis in this paper assumes adiabatic 
behaviour, we have learned elsewhere that universality 
classes of equilibrium systems at their adiabatic transi- 
tions can become universality classes of non-equilibrium 
systems at fast quenches [3D]: these points will be the 
content of a separate paper. Other works [18, 22, 24, 28 
based on studying the dynamics of string networks with 
junctions also suggest that at late times only the lightest 
strings may remain. 



The paper is set up as follows. In section [IT] we first re- 
view some relevant aspects of string statistical mechanics 
in the simplest case: one type of string and no junctions. 
In particular, the duality between strings and fields is 
discussed. In section III we still consider only strings of 



a single tension and type, but now these arc allowed to 
meet at a junction. This section paves the way for sec- 
tion |IV| in which we consider the general case of strings 
of three different tensions o~ a and types, meeting at junc- 
tions. 



As explained in scction [III] there is significant complex- 
ity involved in adding junctions when discussing string 
statistical mechanics, and hence this section is central to 
the development of the paper. Furthermore, technically, 
junctions can be introduced in different ways, and as a re- 
sult we are forced to discuss in detail two specific models 
('bosonic' and 'fermionic') to do so. While bosonic mod- 
els are closer to the physical system we eventually wish 
to describe (and discussed in section [n] when there are 
no junctions), only fermionic models can be generalized 



to the three-string case of section IV At the end of sec- 
tion [TllJ we compare these two models, and conclude that 
they both essentially agree in their phase structure. This 
justifies the use of fermionic models in section [TV] where 
the analysis resulting in the conclusions drawn in figure [JJ 
is straightforward. Finally, we also show, following ideas 
from QCD, that the string system with junctions can be 
rewritten a generalised spin model (XY model). 



II. UNDERSTANDING THE HAGEDORN 
TRANSITION 



In this section we discuss the nature of the Hagedorn 
transition for strings of a single type, with tension a, and 
no Y-junctions. 

As mentioned in the introduction, we proceed by us- 
ing the duality between string configurations and fields 
to write the partition function for the string network as 
that of an effective field theory [31]. As a result, the 
Hagedorn transition can be mapped onto a transition of 
the effective field. Furthermore, provided the right ques- 
tions are asked, one can work with the canonical rather 
than the microcanonical ensemble. 

Consider a classical static picture of non-interacting 
strings in _D-spatial dimensions. These are taken to lie 
on a hypercubic lattice of spacing a, and the energy E 
of the strings only depends on the total string length L 
through E = aL. Near the critical temperature, corre- 
lations are large and the details of the lattice structure 
should be unimportant. We also assume that the net- 
work can be thought of as a set of random walks. Now 
recall the duality between (non-oriented) Brownian paths 
in D spatial dimensions and a scalar field ip of mass to, 
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as exemplified by the identity 

( v (xMo)) = 

f°° 2 /- X ( T )= X 

/ dTeT Tm / Pxexp 

JO Jx(0)=O 



dr'U^ 
o A^' 



This identity can be used to construct an effective action 
(or, more accurately, a free energy) for the string parti- 
tion function Z at temperature T — /3~ x in terms of tp as 

m 

Z = I T>ip exp 



— dx 



D 



4D V r ' 2 



, (1) 



where the mass term is 
M 



The Hagedorn transition temperature, T# = Pjj 1 , is the 
solution to 



J(0) 



-{3<TU 



1 

2D' 



(2) 



The normalisation of ip has been chosen here so that M 2 
is dimensionless. (Note that one would have recovered 
the same temperature Tjj for a gas of strings by counting 
single-loop configurations on the lattice [3J[5]). 

It is important to observe that below the Hagedorn 
transition T < Tjj, ip is a massive free field with M 2 pos- 
itive. For T > Tjj, with M 2 < 0, it describes a tachyon. 
Here fluctuations are large and for this reason the canon- 
ical ensemble often dropped in favour of the microcanon- 
ical ensemble |10j . However, in the conventional picture 
of spontaneous symmetry breaking we arc familiar with 
the way in which tachyons describe instabilities (in field 
space) ; they are understood as corresponding to an inap- 
propriate choice of ground state, the true ground states 
appearing naturally once back-reaction is taken into ac- 
count. 

For example, the inclusion of a repulsive point- 
interaction modifies the free energy to |31j 



S = dx 
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(3) 



thus permitting (tp) to remain finite for T > Tjj. For our 
(p, q) networks, the system has more complicated inter- 
actions than such a simple local repulsion. In particular, 
were the strings allowed to interact at Y-junctions, we 
would expect them to induce additional cubic flip 3 terms 
— as we shall sec in a different context below. However, 
the general implications are much the same. 

The vanishing of the order parameter (ip) at T < Tjj 
can be understood in the following way. Examination of 
the partition function shows that total string density is 
proportional to (<p 2 ), whereas (ip) 2 measures the density 
in infinite string (i.e. string that crosses space) [3"Tl 132) . 



It is the vanishing of infinite string that characterises the 
Hagedorn transition, and not the vanishing of string. 

Although large loops are energetically unfavourable, 
some loops will always exist below the transition (in an 
adiabatic limit). Superficially, free energies like ^ look 
like those of high-temperature quantum field theories on 
dimensional compactification. Either from calculating 
the thermal propagator for excitations at the relevant 
groundstate or by counting microstates of a loop gas we 
get the same result that, in the vicinity of the transi- 
tion, the loop distribution is dominated by the smallest 
possible loops (the ultraviolet limit) 32J. 



III. MEAN FIELD TRANSITIONS; XY MODELS 

As discussed in section [TTJ we anticipate that Y- 
j unctions will induce cubic interaction terms in the dual 
field theory. However, we do not know how to introduce 
them in the exact framework of section |Hj even when 
the junctions are between strings of the same type and 
tension a — the setup considered in the present section. 

In this section we discuss a mean-field procedure which 
allows junctions to be incorporated, and which shows how 
such cubic interaction terms arise. As in section [TTJ one 
can then construct an analogue effective potential, V(tp), 
for a field tp, whose vanishing describes the transition. 
Unfortunately, it is not possible to extend this construc- 
tion to the full effective action, and as a result it is not 
possible to identify the field fluctuations that describe 
finite loops: our analysis is restricted to infinite string 
and the transitions triggered by its creation. Nonethe- 
less, knowing that loops are there enables us to complete 
the picture, qualitatively. It is the mean field procedure 
presented in this section which will be generalised to the 
three-string model in section |IV| 

Again we work in D spatial dimensions, on a periodic 
hyper-cubic lattice of N sites and lattice size a = 1. Let 
i label a lattice site, and /j, = 1, . . . , D the (positive) unit 
vectors in D-dimensional space. There is now a technical 
complication, related to how we allow the strings to pop- 
ulate the lattice. Although there is an energetic penalty 
in having more than one string on a link, in the first in- 
stance we do not wish to restrict the number to unity. To 
do so could imply an effective repulsion between strings 
that is a lattice artefact, and which might induce mis- 
leading terms in the effective potential for the analogue 
field ip. Without this restriction the models are termed 
'bosonic'. 

Models in which at most one string (of any type) can 
lie on a link are termed 'fermionic'. In practice we shall 
find, when we come to mimicking (p, q) strings, that only 
a fermionic model can accommodate junctions of three 
string types. 

An important result of this section is that our concern 
about fermionic models is largely unjustified (though we 
feel it is necessary, for reasons of clarity, to discuss it 
in detail) : both bosonic and fermionic models essentially 
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agree for the small ip values that are relevant for tran- 
sitions, and for which the mean field approximation is 
more reliable. Further, both of these models rewrites the 
string system as a generalised XY model, permitting us 
to think of the Hagedorn transition as one of spin order- 
ing. This suggests ways of going beyond the mean field 
approximation, although we shall not do so here. 



A. Bosonic models 

With conventional lattice notation, let nf be 
the number (0, 1,2,.. .) of strings (anti-strings) on the 
link between the lattice points i and i + fi. 

For strings with no junctions, the Hamiltonian 



A? 



H 



(4) 



gives the requisite energy E — aL to a network of total 
length L. 

Now, depending on the string network we wish to 
model, there is more than one way to proceed. We dis- 
cuss the mean field potenial in each case, making links 
with sections HT1 and ITVl 



1 . Massless junctions 

First we allow the strings to have TV^-fold massless 
junctions i.e. no extra cost in energy. [We are primar- 
ily concerned with N v = 3.] Since the junctions consid- 
ered are massless they do not appear in the Hamiltonian, 
which is still given by Q. 

Rather, the existence of junctions imposes constraints 
on the nf (n~^). Junctions or anti-junctions are permit- 
ted on site i provided the flux into that site is an integer 
multiple of N v : 

a i = £ [( n t» ~ - " n 7-^)] = mod N v, 

(5) 

a constraint which can be implemented through 



1 N " 

S a =o mod N v = /"^ eia9i where B l 



2irki 



Using this representation in the canonical partition func- 
tion 



Z 



cti— mod N v 



enables us to write Z as 



£< 



V 



J 



where the different signs in front of the lattice variables 
9i in the two terms in round brackets reflect the signs in 
([5]). The summations can be performed, to obtain 



Z 



N 



- E,„ ln(l+ J(/3) 2 -2J(/3) cos(9*-9«- ff )) 



where J((3) — e as in (|2J. That is, the Hamiltonian 
of the network is, up to a constant, 

&B = ]Tln[l + J(/3) 2 - 2J(/3)cos(0 i - 6^)}. (6) 



It is not possible to evaluate Z exactly. Hence we resort 
to the mean field approximation scheme (see for example 
|1J), which consists of introducing a trial Hamiltonian Hq 
in which each variable of the system is decoupled from 
the other but depends on an external constant source ip. 
An obvious choice here is 



H {p) 



£< 



(7) 



On writing 



H = H (<p) + [H-H (<p)], 



then 



*= £ 



e -f}H (<p) e -f}[H-H { V )\ 



config 



where the zero subscript denotes (^-dependent averaging 
with regard to H (<p)- As a result the free energy F — 
—Tin Z satisfies 

F{p) < NV{p) = F (p) + (H) - (H ( V )) , (8) 

where V(p) is the mean field effective potential (and F — 
— T In Zq) . Our aim is then to minimize V in order to find 
fmin, which determines the density of infinite string (see 
below) . 

We now carry out the calculation explicitly in the case 
of Y-junctions for which N v = 3. Then 



Z (ip) = 




Q (£> cos(27rfc/3) 



fe = l 



N 



N 
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where 

Now use the results that 

00 m 

(ln(l+p 2 -2pcos6»)) = -2 V —(cos 

z — ' m 



i mb 



(bl < i) 

(9) 



for all measures, and that 



(cos m0) 



(10) 



for the case in point, where 

I m {<p) = i^e vcos(27rfc/3) cos(2^mfc/3) 
3 k 

= \ > {e lp + 2e~ v/2 cos(27rm/3) 



is a discrete version of the Bessel function. Hence, using 
m\i we obtain 



(3V(<p) = -ln(Jo(¥>))+V5 



-273 £ 



m— 1 
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J(f3) m ( 




m I 





(11) 



The periodicity (modulo 3) of the ImW) enables us to 
perform the summation explicitly, to give 



PV(tp) = -]n{e v + 2e-^ 2 ) + cp(^ 
-2DG((3) 



-<p/2 



e <P _ e -(p/2 



where 



1 f l + J + J 2 
G= 3 1 H (1-J)2 



e 1 ^ + 2e-^/ 2 



(12) 



for small J(/3). 

Notice that, because the sum over m in (111 just repro- 
duces the first term with a modified coefficient, V(tp) of 
( 12 1 can be shown to be exactly the mean-field potential 



arising from the Hamiltonian 

GQ3) 



Tjdisc 

H XY - ~- 



E 



(13) 



i.e. the Hamiltonian for a system of unit spins in the plane 
with nearest neighbour interactions in which their rela- 
tive angles are constrained to multiples of 2tt/N v (here 
N v = 3); a discrete XY model. The mean field trial 



Hamiltonian Hq in this case is Hq(p) 



jn. J2 t Si for 



an arbitrary unit vector n in which the spins are decou- 
pled; in other words, an external magnetic field propor- 
tional to (p. 

In order to understand the phase structure of the 
model (either as a spin system or as a gas of strings with 
junctions), consider first the series expansion of V(<p); 



13V (ip) 



1 



2 2 
m ip 



1 



1 



2 r 3 
up to constant terms, where 



flip 



-M 4 



(14) 



m 2 = ^(l-2DG),fi=^(l-3DG),\ = -* 6 (l-2DG). 

(15) 

Observe that the field becomes massless at the tempera- 
ture for which 2DG(f3) = 1, which is in good agreement 
with the Hagedorn temperature of the free dual theory 
of (JIJ since G(j3) ~ J(/3) for J = 1/273 < 1. Further- 
more, as anticipated, the Y-junctions have induced a cu- 
bic term in the potential. In addition they have also 
induced a quartic interaction, vanishing when the field 
becomes massless, that is repulsive when the field be- 
comes tachyonic. 

As a result of the cubic term, the potential in equa- 



tion ( 12 1 can be shown to have a weak first order phase 
transition. The critical temperature, however, cannot be 
obtained from (15) as it occurs at values of ip ~ 1. Nu- 
merically, however, one finds that 2G cr i t (D = 3) ~ 0.31 
and 2G cr it(D — 4) ~ 0.23. We shall not consider the 
first order transition further, since it is not reliably ro- 
bust against rapid quenches which is what we ultimately 
have in mind. 



2. Massive junctions 

Alternatively one might want to model string networks 
with massive junctions — that is, is to introduce junc- 
tions with an energy cost v. (These can model massive 
monopoles, which may be formed at the vertex in dif- 
ferent symmetry breaking schemes [33 .) We can then 
recover massless vertices by taking v — > 0. Furthermore 
this construction allows one to calculate the average den- 
sity of vertices at temperature T, by simply differentiat- 
ing Z with respect to v. This will be discussed in section 

El 

To add massive vertices, we allocate a vertex number 
(0, 1,2...) to each lattice site, constrained by 

a i = E [( n t,M ~ n t-w) - K> - n i-w)] 

+ 3{pt- P r)=0 (16) 



for Y-junctions, while the Hamiltonian acquires an extra 
term 



P 



± 



N 



(17) 
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Performing the sums over the nf and the pf leads to a 
Hamiltonian 

m = -^ln[l + J 2 (/3)- 2 J({3) cos(0 i+/1 - 0;)] 
- ^ln[l + A 2 (/?) -2A(/3)cos36y, (18) 

where the Oi are now continuous variables, the Lagrange 
multipliers that arise from imposing the constraints 



fin^n — 



1 

2tt 



Also we have defined 
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(19) 



(20) 



analogously to J in Then, carrying out the same 

mean field treatment as above yields 



0VW{<p) = -ln(/ (¥>)) + ¥> 



Io(<p) 



-2D £ 

m=l 

" K(f3) m fl 3m (<p) 



-2E 

m— 1 



(21) 



where the /,„ are (continuous) Bessel functions. 

For non-zero A cubic terms arise from the A3 Bessel 



function, to give rise to a potential of the form ( 14 1, with 
coefficients 



1 



(1-2£>J), fi=- 



A 



A = — 



1G 



(22) 

As expected, we have tachyonic instability at J = 1/2D 
and a cubic term in the potential. 



The slightly different behaviour of (22 1 and (15) is to 



be expected, since since we arc implementing the bound- 
ary conditions that count vertices differently in the two 
cases: in other words, they correspond to different imple- 
mentations of the mean field approach. However, since 
the mean field result is, strictly, an upper bound, we 
could, if we wished, only retain that solution that is nu- 
merically lower. In practice, this is not necessary since 
there is close numerical agreement at relevant tempera- 
tures. Massless junctions correspond to taking A = 1 for 
which /i = — 1/8 , the value arising in (151 when 2DG = 1 . 



Further, a numerical study of (21 1 shows that the transi- 
tion tends to become first order as A — > 1, in agreement 



with the discussion of (15 1 



3. No junctions 

For continuous 'bosonic' string with no junctions both 
approaches give the identical result. In the first case, 



we eliminate junctions by taking N v — > 00, whereby the 
discrete Bessel functions are replaced by their continuous 
counterparts. In the second, taking v — > 00 (A = 0) just 
recreates the same series. 

In each case, on expanding V^°\(p) for small tp we 
recover the second order transition at the Hagedorn 
temperature Tjj of Section II (see equation (|2])) when 
2D J ((3) = 1, and when the ip field becomes tachyonic. 
However, it can be seen that of (21 1 becomes un- 
bounded below as T — > 00. This is not quite the be- 
haviour of ([lj , for which the potential is unbounded be- 
low for all T > Th, showing the limitations of the mean 
field approach for very large \<p\. Nonetheless, this simple 
example shows how the introduction of vertices induces 
interaction terms in the effective potential to stabilise the 
ground states. 



B. Fermionic models 

We now consider the most simple 'fermionic' models. 
It is these which can straightforwardly be extended to 
the general three string-type model of section |IV| We 
will also address the concern raised at the beginning of 
this section: that the 'fermionic' model might add an 
effective repulsion between strings, which could induce 
misleading terms in the effective potential. We will show 
that this is not the case. 

Thus, we now restrict the number of strings on each 
link to ni yfl € {0,±1}. That is, the link from site i to 
i + fi contains either a single string, a single anti-string, 
or no string at all. 



1. No junctions 
With no junctions, the Hamiltonian is 

N D 

tf = EE™L> ( 23 ) 



= 1 jU.= l 



subject to the constraint 



(24) 



Performing the sums over the n i fl leads to a Hamiltonian 
PH = - YM l + 2J ^) cos (^+m - *<)]. (25) 



1 .pi 



where the Oi are again the Lagrange multipliers that arise 
from imposing the constraints 



&cti,0 — 



1 

27 
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d9,. 



Defining J by 



J 



J 



1 + J 2 ' 



(26) 



(27) 
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whereby J {(3) ~ J {(3) when J -c 1, a similar calculation 
to that above (see also section IV I shows that the mean- 
field potential is, for J < 1, 



PV^iv) = -lnlo(p) 



h{v) 
h{<p) 



rn — 1 



(28) 



to the mean-field approximation. In practice, we know 
of no work on the generalised XY models appropriate to 
the three-string models (with or without junctions) and 
stay with the mean-field approximation. 

To give a meaning to (p we note that the average den- 
sity of (infinite) strings is proportional to tp 2 , as antici- 
pated, given by 



The tp field now becomes massless at 2D J {(3) — 1, with 
a second order transition. With J ss J this is slightly 
displaced from that of the bosonic strings but, at the 
qualitative level at which we are working, can be said to 
agree. Note that both potentials (28 1 and (21 1 show a 
Z 2 symmetry under tp — > —tp that is broken above Tjj, 
and restored below Th, contrary to the usual pattern of 
symmetry breaking, but as in section |TT| 

On comparing (28) with (21 ) we see that they differ in 



that the former has alternating signs in the Bessel func- 
tion expansion, whereas the latter does not. Because 
higher terms in the series in powers of J{(3) become sig- 
nificant only at increasingly large tp, the artificial repul- 
sion induced by the 'fermionic' assumption (that is, of 
no more than one string per link) is a large-tp effect in 
the mean field approximation, and hence where the ap- 
proximation is at its least reliable. However, since the 
transitions are determined by small tp, we can use either. 
This is an important result of this section. 

In fact, for J small, both approximate the mean-field 
potential of the XY-model, with spin-spin Hamiltonian 



H 



XY 





23 



%) 



£* 



i+fJ.- 



This follows from expanding (|25j), for which 



V X y{v) = -hi/ (<ys) + ^ 



h{<p) 



-2D J {13) 



(29) 



(30) 
1. Rather 



showing a second order transition at 2DJ{[3) 
than just perform a series expansion in tp as in ( 15 1 , more 



generally we see that extrema of Vxy{<p) satisfy 

tp-ADJ{f3)u{tp) = 0, 
where u{ip) = h{tp) / I {tp) . 



(31) 



tp = is always a solution to (31). For 2D J {(3) > 1 



there is a further pair of solutions, itp, tp > 0, which are 
the minima. We note, for future use, when we need to 



count extrema, that (31 1 behaves like the cubic equation 
obtained from just retaining terms up to 0{tp 4 ) in the 
expansion of the potential in the existence of three roots. 
The inclusion of higher terms in the scries in J does not 
seem to affect this empirically and it is not necessary to 
go beyond the XY model, now and hereafter. 

When the XY model is a good approximation we could, 
in principle, use known results about it without resorting 



1 

w 



-J{(3)(3 



dV- 



XY 



dJ 



ADJ{(3) ' 



(32) 



2. Massive junctions 

We end this section by including Y-junctions in the 
fermionic model (still of a single string type). Given 
that the occupation numbers are limited to 0,±1, there 
is no analogue of the mod 3 description for massless ver- 
tices discussed in the bosonic case (see equation ^ ) . We 
therefore consider massive vertices. There is now a single 
vertex number pi = {0,±1} constrained by 

on = ^ K> ~ n i-fi,fi] + 3Pi = ( 33 ) 



with the Hamiltonian acquiring an additional term 



On defining K by 



K = 



K 



1 + K 21 

the mean-field potential is, for ( J, K < 1), 



(34) 



(35) 



f3vj K \tp)=0V« 3) {tp)+2j2 



{-K{P)) m (hrrM 

m \ I {tp) 



_ (36) 

[This follows from the generalisation of Q, used earlier 
in (28 1 that, up to a constant, 



{\n{l + 2K cos a)) = -2 ^ 



(-KY 



(cos ma) (37) 



for all measures and K < 1, together with the specific 
result 



(cos m6% 



_ J ffe ycose cosme = I m (<p) 

h{f) 



cos 6 



for our choice of Ho.] We note that unfortunately, for 
a simple cubic lattice, the requirement that K < 1, 
necessary for convergence of the series in ( 37 1 , imposes 



K < 1/2. Hence that the mean field approximation is not 
valid for light vertices in the fermionic case (as opposed 
to the bosonic one in (21 1). We consider this constraint 



to be an artefact of the lattice fermionic approximation. 
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Despite that, note that mean field potential (36 1 leads effective potential will take the generic form 



to an XY model in the presence of an external source 
[3J [5] in which we retain only the first term in the power 
series in A" in (36) (or the first term in the series in K in 
([21)). As a result, there is always a second-order transi- 
tion, as in the bosonic case. 



Finally we also note that the density of string (32 1 is 
unchanged by the inclusion of junctions. 



C. Summary of section |III| 

In summary, in this section we have seen how the inclu- 
sion of Y-junctions in a model of a single string type can 
provide the back-reaction necessary to prevent tachyonic 
instability at the Hagedorn temperature. Further, pro- 
vided we restrict ourselves just to infinite string, whose 
density is the order parameter, we can go beyond the 
Hagedorn temperature, still with the canonical ensem- 
ble. 

We have also discussed two models, 'bosonic' and 
'fermionic', and shown that the concern raised about 
fermionic models at the beginning of this section is un- 
justified: both models essentially agree for the small ip 
values that are relevant for transitions, and for which the 
mean field approximation is more reliable. 

We have also shown how the value (p of the field at 
the minimum of the effective potential is related to the 
density of infinite strings in the system. As we discuss 
in section |IV| it is equally apparent that the density of 
vertices is also determined by (p and obtained by differ- 
entiating the partition function with respect to v. 

With this behind us, we now consider the case of three 
different string types with Y-junctions, as a model for 
(p, q) strings. We note that, oddly, the analysis of QCD 
confinement of [3J HI [S] , that we have called upon in this 
paper, was performed in the context of a single-string 
model, not permitting 'colour'. Although this was not 
our intention, a more realistic description of QCD is given 
by the model that follows, in the limit of equal tensions, 
in which our three string types correspond to coloured 
flux tubes. 



IV. THREE STRINGS, FERMIONIC MODEL 

The basics of our model are the following. 

As stated in the introduction, we model the (p, q) 
string network by a network of three different types of 
fundamental strings, labelled by a = 1, 2, 3 as red, green 
and blue, say. Generally the strings also have different 
tensions a a . The strings do not interact with each other 
(nor with themselves), except at a Y-junction (or vertex) 
which is defined to be a point at which three strings of 
different colours meet. 



(3V(<pi,<p2,<pa) = E 



2 1 



1 



2"" ra ' 4 
+ H(piip 2 f3 



\ 4 



(38) 



Potentials of the type ( 38 1 , with temperature-dependent 
coefficients, have been studied in other contexts 
e.g. transformations of vortex types in superfluid 3 He 



We know that (38 1 is valid if Y-junctions are excluded, 



when fi = 0. In this case, from the single string models 
ml<x(l-2DJ a (fi)) ) (39) 

with J a — e~ aal3 . In the following discussion we suppose 
that 



01 < C2 — CT 3 

The critical Jg rit = 
temperatures /3 a 



Ji > J 2 > J3 



(40) 



= 1/(21?) define three critical inverse 
T- 1 with 



03 < fa < Pi 



(41) 



in the vicinity of which mr a cx (1—T/T a ). That is, with no 
interactions we expect three sequential Hagedorn transi- 
tions as, on cooling, the heavier strings disappear from 
the picture, leaving the lightest until last before it disap- 
pears in turn, leaving just small loops. 

Our aim is to understand the effect that Y-junctions 
have on this picture. 



In practice, we are not able to recreate ( 38 1 in a bosonic 



model with coloured Y-junctions, with arbitrary numbers 
of strings on each link. (The reason is that we are unable 
to write down a generalised form of the constraint ([5| in 
the 3-string case.) We therefore restrict ourselves to a 
fermionic model, in which there is at most one string of 
each type on a link. As discussed in the previous section, 
we expect that the effective repulsion this implies can be 
ignored at small field values. As in the case of the single 
string type, in order to be able to use mean field theory 
we are obliged to give the vertex a non-zero mass v. 

As before, we assume that the energy of the different 
strings is proportional to their length L (E a = cr a L). 
The different strings are described respectively by the 
variables nf , which all take values in {0,±1}. There 
are also vertices, described by the variable pi £ {0,±1}, 
joining strings of 3 different types. The Hamiltonian of 
the system takes the same form as the for the single string 
case, 



* = E 



(42) 



We now need to impose the constraint that a junction 
is where three different colour strings meet: this is done 

by 



Following equation (14 1, our expectation is that the 



7; 



r = E«a 



1— /J.,/X 



^)+ Pt = 0, Va. (43) 
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Although summing over a would essentially recreate the 



constraints (33 1, equation (43 1 is more specific. In par- 



ticular, (43 1 does not forbid different string types from 



lying on top of each other. 

As in the previous section, the constraints are imposed 
in the standard way through Lagrange multipliers, which 
is equivalent to writing the Kroneker delta as 



#Tf,o — 



1 

2^ 



271 



d0f 



(44) 



(no a summation) for each 7" . Hence the partition func- 
tion is 



Z((3,v,a a ) 
dO? 



i, a 



x > e 

P: 



which, on carrying out the summations gives 
d9 



(45) 



i,a f-L 

x[l + 2Kcos(]TO- (46) 

a 

This corresponds to the Hamiltonian 

0H=-J2 M 1 + 2 J a cos(0f +M - Of)] 

i a 

« £ 2 J a cos(^ +M - 0?) + 2 2* cos(]T 0?) (47) 

for small J a and i^T . 

The mean field treatment therefore contains three vari- 
ational parameters ip a . Following the same steps as in 



section III the trial partition functions which decouple 
different lattice sites are 



Zo(/3,a a ,cp a ) 



n 



d9f 
2tT' 



= [Jo (<?«)] , 
(48) 



while the mean field effective potential is 



a 



- lnl (ip a ) + 



» (-K)™ ( Info) I m (tp 3 y 

^ m ^ /ofai) ^0(^2) Jo(^) ' "" 



where each J a is defined as in (27), and K is given in 
pi. 



As discussed in section |III| it is sufficient for our pur- 
poses to approximate (3V(ip a ) by the first term in the 
scries of (|49|), 



/3V XY (<p a ) = J2 

oc 



- lnl (<p a ) + 



2DJ 



2K 



h((fii) h(<P2) h(<P3) 



Jo(ipi) Io{v>2) hivz) ) ' ^ 
This corresponds to making the small </, K approxima- 
tion in ( 47 1 . That is, the model ( 46 1 is a generalised XY 



model, consisting of three spin-like variables defined on 
each lattice site i, making angles Of with respect to some 
fixed axis, interacting amongst themselves through the 
if-dependent term. 

We have achieved our goal in that, if we expand 
Vxy(<P*) of ((50J (or, indeed the full V{ip a ) of |49])) in 
powers of ip a we recover the generic potential ( 38 1 as the 
first few terms in the series. 

However, we can say more. As in our earlier examples, 
attaching a nominal energy to each vertex allows us to 
calculate the density of vertices. Specifically, the density 
of vertices on infinite strings is 



' dK 

ex (pi<p 2 (p 3 (51) 

at the minimum (^1,^2,^3) of VxY^a)- The small 
loops corresponding to the field fluctuations that are in- 
visible to our mean field analysis contain vertices not 
counted in (51 1. 
As in section J 



III we now look for the extrema of the po- 



tential in order to determine the density of infinite string 
and the density of vertices. As expected from section [Hi] 
a full numerical analysis (that we have performed) with- 
out the XY approximation does not alter our qualitative 
conclusions and barely changes our quantitative results. 

As we noted earlier, in the main works on QCD ([2111]) 
all flux strings were taken to be of a single kind, leading 
to a very different potential, in which Ji(¥>i)Ji(y2)Ji(y3) 
is replaced by I^ip) for example. In particular, as we 
shall see later for ( 49 1 , with equal tensions there is no 



first-order transition when there are three string types. 

A. K = 0: no vertices and three independent spins 

We have already anticipated the results for this simple 
case, but it is helpful to see them in greater detail. For 
K = the XY model reduces to three independent, un- 
coupled, XY models with Z 2 x Z 2 x Z 2 symmetry under 
ip a — > — ip a . The extremal points are when 

dV M F 



dip a 







4DJ a 



u{(p a )=Q (52) 
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where u(ip) = Ii(ip) / Io(ip) as before. One possible solu- 
tion is always (p a = 0, the only real solution if 2DJ a (f3) < 
1. 

If 2DJ a (/3) > 1 then there are two further real solu- 
tions, denoted ±tf a , where we take (p a > 0. The 3 3 = 27 
possible extrema ip — {(p\, (f>2, fz) then break down into a 
non-degenerate <p = (0, 0, 0), three doubly degenerate so- 
lutions, exemplified by ip — (±<pi, 0, 0), three fourfold de- 
generate solutions, exemplified by (±y>i, ±<f>2, 0) and an 
eightfold degenerate solution (±<pi, ±02, ±03)- It is suf- 
ficient to restrict ourselves to the positive sector f a > 0. 

To determine which of these are maxima, which min- 
ima, and which saddle points we need to calculate the 
eigenvalues of the Hessian M 7 s = 9 2 Vxy I dtp^dips at the 
extrema. An extremum is a minimum if all are positive, 
and a maximum if all are negative. Otherwise one is 
dealing with saddle points. 

With K = 0, the only non-zero entries are on the di- 
agonal with (no summation) 



u'(f a )[l-ADJ a (P)u'(f a )}. 



(53) 



For the case in hand the answer is very simple and very 
obvious. 

1. (3 > Pi(> /3a, /3s). lu this range the global mini- 
mum occurs at <p = (0,0,0). 

2. /?2 < P < Pi- Now (0i, 0, 0) is the global minimum. 
[(0,0,0) is a now saddle point.] 

3. /?3 < P < /?2- In this range it is easy to see that 
(0i, 02, 0) is the global minimum. 

4. P < P 3 . Here it is equally straightforward to see 
that (0i,02?03) is the global minimum, (0,0,0) is 
a maximum, and all other points are saddle points. 

As expected, as the temperature is increased infinite 
strings of the lightest tension first are nucleated at P = 
pi; then those of the next lightest tension at P — P2; 
and finally the heaviest strings when P = P 3 . When 
one decreases the temperature from a very high one, the 
opposite happens. 



B. K 7^ 0: vertices and three coupled spins 



Let us now consider the effect of Y-junctions in the 



generalised XY model of ( 50 ) . For unequal a a the sym- 
metry of Vxy is now explicitly broken from Z2 x Z2 x Z2 
to D 2 

Pi 
Pi 
P 3 



'2 x Z 2 , generated by 

<Pi, 



<P2, 



<P3 -i 
<^3 
<^3 



-<P3 
<f>3 



<Pl -> ~<Pl, V2 
<Pl -> ~<Pl, V2 -<P2, 

If any tensions are equal the symmetry is correspondingly 
increased. Imposing 8Vxy /d(p a — gives (no summa- 
tion) 

u'^ a )[^ a -ADJ a (P)u^ a )-2K(p)u(ip p )u(^)} = 

(54) 



where P — (a + l)mod3, 7 = (a + 2) mod 3. There 
are obvious solutions to these coupled equations: (0, 0, 0) 
for all P; (fi,0, 0) with ipi — (p\ (the standard solution 
provided 2DJ\ > 1). The important point though is that 
it is not possible to have a solution with only, say (p~\ = 0, 
and the other two non-zero. One can see this from (54 1, 
where setting (pi = would require that one of the other 
two y>'s must vanish. 

At the extrema the Hessian has the same diagonal el- 



ements as in ( 53 1 , but off-diagonal elements 



M ajj = -2K(P)u'((p a )u'(ipf3)u(ip 7 ) 



(55) 



We now evaluate these at the different extrema identified 
above and discuss the consequences. 

Case 1: (p a = 0,Va. 

This reduces to the free-string case above, as here 
the off diagonal terms of M also vanish. We have a 
global minimum for P > Pi as all the eigenvalues are 
positive. Otherwise, when @ 3 < P < Pi we have a saddle 
point, and for P < P 3 a global maximum. 

Thus, as the temperature increases (or P decreases) 
the 1 direction will 'roll' first. 

Case 2: (p2 = ^3 = but (pi ^ 0. 

Now, notice that the temperatures j3±, P 2 and ^3, 
as defined for free strings, are in principle relevant only 
when ip a = since then the off-diagonal terms of M 
vanish. When non-zero ip a enter, we have to worry about 
the off-diagonal terms, and find the new eigenvalues. 
This in turn will introduce new critical (if-dependent) 
temperatures. 

As before, (pi is the solution of the standard equation 
provided 2DJi > 1 or /3 < P\. 

When P — P 2 the smallest eigenvalue is negative, show- 
ing that (<pi,0,0) is not a local minimum. There is an 
intermediate temperature the solution to 

(1 - 2DJ 2 (/3*))(1 - 2DJ 3 (P*)) = K 2 {P»)u 2 {fi) (56) 

that denotes the transition from local minimum to saddle 
point. That is, strings of type 2 and 3 are nucleated at 
the same time. 

To summarize: for P > Pi there is a global minimum 
at f a = 0. For p 2 < P* < P < Pi the global minimum is 
at (01,0,0). 

Case 3: (pi,if2,(p3 all non-zero. 

For type 2 and 3 strings are nucleated 

since one cannot have only one non-zero ip a . Hence we 
expect to have non-zero if a for all a. However, there is 
nothing at this stage to preclude the possibility of even 
further transitions, of first and second order. 

Discussion. 
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cp 2 



FIG. 2: Schematic representation of the trajectory of (p in 
field space. The arrow indicates the trajectory as a function 
of decreasing temperature. 



it stands the analysis above is restricted to closed or 
infinite string. The addition of quarks to string ends 
changes the picture again. Further, since flux tubes are 
not fundamental in any sense, the 'Hagedorn' transition 
in QCD has a different status, with no ambiguity about 
increasing the temperature beyond it. 

Density of vertices: 

Finally we end this section with a comment on the 
density of vertices in the different phases. From (51 1, 
and since 7 m (0) = for m > 1 it follows that, on 
differentiating V(ip a ) with respect to K, 







(57) 



when any (p a = 0. Thus, we only have a non-zero 
density of vertices on infinite strings for f3 < /3*, i.e. at 
temperatures high enough for all infinite string of all 
types to be present. This is shown in Figure [l] 



We can get some help from elementary Morse the- 
ory, applied to the whole (p a space and not just the 
positive sector [34 . Empirically, for the purpose of 
counting extrema, equations (54) also behave just like 



the cubic equations that would follow from taking only 
the leading terms (38 1. According to this, when we 
have 27 extrema, no more than 14 can be minima. The 
cases of all (p a = or one Cp a non-zero may produce 
7 (7 = 1 + 3x2) real extrema and therefore 20 may 
correspond to extrema with no (p a vanishing. From D2, 
each is fourfold degenerate, implying that there may 
exist five (5 = 20/4) different least symmetric extrema, 
of which no more than three can be local minima. This 
still allows for either first or second-order Hagedorn 
transitions as (3 is reduced below /3* (or temperature 
increased) . 

Now consider the case when two string types have (ap- 
proximately) the same tension, and the other is markedly 
different, e.g. one string is very light, and the others 
heavy. The cases of all <p a = or one <p a non-zero still 
may produce 7 (7 = 1 + 2 + 4) real extrema. However, 
each extremum with no Cp a vanishing is now approxi- 
mately eightfold symmetric. As a result we do not expect 
more than two of them, of which only one can be a local 
minimum. This means that there cannot be any further 
transitions as f3 is reduced below /?„ . Although a first or- 
der transition cannot be precluded, empirically we have 
only found second order transitions even for a a taking 
different values. The situation is summarised schemati- 
cally in figure [2] 

From the above discussion, it follows trivially that, 
for equal a a , (with twelve-fold degeneracy for all (p a 
non-zero) there is just one second-order transition. This 
is relevant to an idealised version of QCD. However, as 



V. CONCLUSIONS 

The main idea of this paper has been very simple: 
that we can describe the thermodynamics of a network 
of strings of three different types (and tensions) by an 
effective three-field theory whose potential V((px,<P2, ¥3) 
takes the form 



1 2 2 . 1 \ 4 

-m a (p a + -\ a tp c 



A*</?1^2</?3 



(58) 



The interaction coefficient /z reflects the presence of Y- 
junctions at which one string of each type meet. The 
coefficients are temperature dependent, with rr? a oc (1 — 
T/T a ) in the vicinity of its zero. If fx were zero, the 
T a would be Hagedorn temperatures for the individual 
string types. As a result, the discrete symmetries of V are 
broken at high temperature, restored at low temperature, 
in a reversal of the usual pattern. 

Our main results, summarised Figs. 1 and 2, essentially 
follow from the form of ( [58] ) alone, supplemented by an 
understanding of the order parameters, that they char- 
acterise infinite string, and not loops. In consequence, in 
a network of strings of different tensions it is the lightest 
strings whose infinite strings survive last after Hagedorn 
transitions, and even those disappear in turn, to leave a 
collection of small loops. This is despite the presence of 
junctions between strings of different types. That is, the 
only role that the junctions play is in these small loops 
of string whose presence is the only memory of the initial 
proliferation of strings of all types. 

The burden of this paper has been to provide a model 
in which we can see how the potential ( 58 1 is realised, 



almost as proof of principle. This has turned out to be a 
non-trivial task and the model at hand, an extension of 
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similar models used in QCD in a much more restricted 
situation, has its faults. A well as picking a path through 
the 'fcrmionic' lattice artefacts, as in the calculations for 
QCD strings, our strings are also assumed to be non- 
interacting and static. Furthermore we are often pushed 
to consider the model in a limit of parameter space where 
approximations are not always well controlled (just as in 
[3J HI [5]). Our one string bosonic model demonstrated 
how, for a single field, /j,(p 3 terms arise naturally. How- 
ever, being unable to generalise the bosonic model to 
three string types, we have also had to introduce massive 
vertices in the three string model as an artefact of the 
lattice mean-field approximation. Naturally, any specific 
model will give more information than just the leading 
terms of V of ( 58 1. In our case the model is a generalised 



XY model, in which transitions are seen in the language 
of spin ordering and which, in principle, permit better 
than the mean-field approximation. 

As suggested above, our analysis points to the the fi- 
nal stage of the transitions as being that of a single string 
type, collapsing into loops, which was the original case 
to be studied, primarily in the context of Nambu-Goto 
strings. In that case, the full statistical mechanics has 
been studied in detail, and can be generalised to non- 
static strings. The result however, is the same! Indeed, 



rather than consider random walks in space, one can 
consider simultaneous independent random walks on the 
Kibble- Turok spheres for left and right-moving modes re- 
spectively [12]. The microstate density at the transition 
is the square of that for simple random walks, but inte- 
grating over centre-of-mass coordinates reduces the state 
density to that of (appropriately defined) single static 
random walks. 

Another way to make this adiabatic picture dynamical 
is to attempt to determine the time scales of the string 
network transitions from the timescales of the effective 
field theory, using the Kibble scenario [3D]. This relies 
on little more than causal bounds, and the analysis is 
under way. 
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